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ABSTRACT 



In this paper we present the detection of a 1.8 Mj up planet in a 538 day period trajectory orbiting the metal-poor star HD 171028 
([Fe/H] =— 0.49). This planet is the first to be discovered in the context of a HARPS program searching for planets around metal- 
poor stars. Interestingly, HD 171028 is one of the least metal-poor stars in the sample. This discovery is placed in the context 
of the models of planet formation and evolution. 
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1. Introduction 

The discovery of more than 220 extra-solar planets or- 
biting solar-type start0 is providing crucial evidence for 
the processes of planet formation a nd evolution (for a 
review see e.g. lUdrv fe Santosll2007lK Among these, im- 
portant clues come from the study of planet-host stars. 
The well known strong correlation betwe en the presence 
of giant planets and the stellar metallici ty dGonzalezlll99 
Gonzalez et~aT1l200ll:[Santos et al.l l200lll2004bll200alRei 
2002t iFischer fc Valentil 120051 ) suggests that giant plan 



of detected transiting planets orbitin g bright stars (e.g. 
Sato et alJEooi iBouchv et alJl2005bh . 



Interestingly, however, several giant planets were found 
to orbit metal-poor sta rs, with down to ~2 times less met 



als than the Sun (e.g. ISetiawan et al. 20031 : Mayor et al. 
2004t ICochran et alJl20oC Some programs to search for 



planets around such objects were also started. Two of 
these include the use of both Keck and HET telescopes 
([Sozzetti et al.ll2006t ICochran et al.ll2007l respectively). A 
third one, presented in this p aper, is part of the HARPS 



ets are more easily formed around higher metal-content 
stars . This lends supports to the co re-accretion model 
(e.g. lMizunolll980t|Pollack et al]|l996l) as the main planet 
form ation mechanism (e.g. Ida fc Lin 2004ri Benz et al 



GTO planet search program (jMavor et al. I l2003h . 



l2006h . in opposit i on to the alternative d isk-instability 
model (|Bosslll997l 12003 iMaver et al.ll2002t) . 



The higher probability of finding giant planets orbit- 
ing metal-rich stars prompted a number of different sur- 
veys dedicated to metal-rich samples (e.g. iTinnev et al 



2002 1 : IFischer et al.ll2005l ; lda Silva et~alll2006bl : iMelo et al 



20071 ). Given their observing strategy, these programs un- 



veiled mostly short period planets, strongly biasing the 
known samples, while positively increasing the number 



Part of the goal of these programs is to try to under- 
stand how frequent are giant planets orbiting metal-poor 
stars, and what is the metallicity limit below which no gi- 
ant planets can be observed. Such constraints would give 
important new clues about the processes of p lanet for- 
mation and evolution (e.g. iMatsuo et al. 2007 ). The re- 
cent finding that the metallicity correlation may no longer 
be valid for Neptune-mass plane ts (see case of HD 4308 
([Fe/H]=-0.31) and discussion in llJdrv et all (|2006l )). to- 
gether with theoretical predictions suggesting that very 
l ow mass planets may be common ar ound metal-poor stars 



(llda fc Linll2004at iBenz et "all l2006h . renewed the interest 



* Based on observations collected at the La Silla Parana 
Observatory, ESO (Chile) with the HARPS spectrograph at 
the 3.6m telescope, under the GTO program 072.C-0488. 



See e.g. table at http://www.exoplanets.eu/ 



for these surveys. 

In this paper we present the first detection of a giant 
planet orbiting one star from the HARPS metal-poor stars 
survey. The planet orbits HD 171028 ([Fe/H]=-0.49) in a 
538 day period orbit. In Sect. [2] we describe the sample 
used. In Sect[3]we present the observations of HD 171028, 
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Fig. 1. Left: Metallicity distribution for the whole HARPS 
metal-poor sample (photometrically derived metallicities). 
Right: M„ vs. b — y diagram for the stars in the sample. 



providing the stellar parameters and the orbital solution 
found. We conclude in Sect. [4] 



2. The HARPS metal-poor sample 

The HARPS GTO program started in October 2003 
to follow several different samples of solar-type stars 
( Mayor et al. I l2003h . The remarkable long term precision 
of HARPS allow ed the discovery o f several planets am ong 



the targets (e.g. iPepe et al.lbool lLovis et all [20051. 



eluding the large majority of the known planets with 
masses of the orde r of the mass of Neptune or below 
dSantos et al.ll2004aHBonfils et al.ll2005ll2007t lLovis et aL 
l2006HUdrv et alJl2006L l2007h . 

To explore the low metallicity tail of the planet-host 
stars distribution, one of the samples currently studied 
with HARPS is constituted of 105 metal-poor or mild 
metal-poor solar-type stars. This sample was chosen base d 



on the large FGK-catalogue of INordstrom et al.l (|2004l ) 



From this catalogue, we took all late-F, G, and K stars 
(b — y >0.330) south of +10° of declination and having a 
visual V magnitude brighter than 12. From these we then 
excluded all known visual and spectroscopic binaries, all 
stars suspected to be giants, and all those with measured 
projected rotational velocity v sini above ~6.0kms _1 (to 
indirectly exclude the most active stars). Finally, we con- 
sidered only those targets with photometric [Fe/H] be- 
tween —0.5 and —1.5. 

In Fig.Q]we plot the metallicity distribution of the sam- 
ple, as well as a M„ vs. b — y diagram. The distribution 
is clearly peeked between —0.5 and —0.8, while a long 
"flat" metallicity tail is present down to [Fe/H] =—1.5. We 
note that HD 171028 (spectroscopic metallicity of —0.49), 
was taken into the catalogue only because its photometric 
metallicity is lower than this value (see discussion in the 
next section). 

The final 105 stars in the sample have their V magni- 
tudes between 5.9 and 10.9, distributed around an average 
value of 8.7. After a 15 minutes exposure these magni- 
tudes allow to obtain a S/N high enough to derive radial- 
velocities with a precision better than 1ms -1 for the ma- 
jority of the targets. 



3. A new planet around HD 171028 

3.1. Observations 

We obtained a total of 19 measurements of HD 171028 in 
GTO time (program ID 072.C-0488) using the HARPS 
spectrograph (3.6-m ESO telescope, La Silla, Chile). The 
observations were carried out between October 2004 and 
April 2007, and the radial-velocities were obtained using 
the latest version of the HARPS pipeline. We refer to 
IPepe et aL (2004) for details on the data reduction. 

Exposure times varied between 260 and 900 s, and the 
individual photon-noise error in the radial-velocities was 
always below 1.7ms (median of 0.8ms" 1 ). In this error 
we do not quanti fy the uncertainties d ue to the stellar 
oscillation modes ( Bouchv et al.l l2005a[ ) . This "noise" is 
typically averaged out if a long (15 minutes) exposure is 
done. Unfortunately, this was not always the case when 
observing HD 171028 (only 12 out of the 19 measurements 
were done following this strategy). 

The individual spectra were also used to derive both 
the Bisector Inverse Slope (BIS) of the H ARPS Cross 
Corre lation Function (CCF), as defined bv iQueloz et al 
(|2000l ). as well as a measurement of the chromospheric 
ac tivity index logi?^ K , following a similar recipe as used 
bv lSantos et~aU ^2000) for CORALIE spectra. Finally, the 



combined high S/N HARPS spectra were analyzed to de- 
rive stellar atmospheric para meters and iro n abund ances 
using the method described in lSantos et al. ( 2004rJ ). 



3.2. Stellar characteristics 

There is little information available on HD 171028 
(BD+06 3833, V=8.31, B- V=0.61), and only two papers 
appe ar in a Simbad query (IOlsenl(l994l : INordstrom et alJ 
2004). The star is also not in the Hipparcos catalogue 



(lESAlll99Tt> . and most of the information available on the 
Tycho catalogue is not very accurate (e.g. the parallax 
listed is 9.1±7.8milliarcsec). 

The analysis of our combined HARPS spectra, with 
a total S/N ratio above 500, provide T e // = 5663 ± 20, 
log 5 = 3.84 ±0.05, and [Fe/H] = -0.49 ± 0.02. These 
values were obtained following th e methodology and line- 
lists used in ISantos et al.1 ( 2004bh . Similar values for the 
stellar parameters are ob tained using the la rger Fei and 
Fell line-list presented in lSousa et al.1 (|2007| ). and the au- 
tomatic ARE93 code for line-equivalent width measure- 
ment (T eff = 5693±16, log g = 3.85±0.05, and [Fe/H] = 
— 0.48±0.01). For the rest of the paper we will consider the 
first set of parameters. The errors mentioned above denote 
internal errors only. Systematic errors affecting e.g. the 
temperat ure scale are thoroughly discussed in the litera- 
ture (e. g.lSantos et al.ll2004bl (Ramirez fc Melendezl 12004 



Casagrande et all 20061 iMasana et all 2006h . and no 



sensus seems to exist at this point. This discussion is, how- 
ever, out of the scope of the current paper. 



http: / / www.astro.up.pt /^sousasag/ares 
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Table 1. Stellar parameters for HD 171028. 



Parameter 


Value 


Reference 


Spectral type 


GO 


Simbad 


m v 


8.31 


Simbad 


B-V 


0.61 


Simbad 


b-y 


0.43 


Nordstrom et al. (2004) 


log R'hk 


-4.92+ 


This paper 


Distance [pc] 


90 


This paper 


v sini [km s _1 ] 


2.3ft 


This paper 


T cff [K] 


5663 ± 20 


This paper 


logs 


3.84 ± 0.05 


This paper 


€t 


1.32 ± 0.05 


This paper 


[Fe/H] 


-0.49 ± 0.02 


This paper 


Mass [M Q ] 


0.99±0.08 


This paper 



f From HAR PS spectra u s ing a calibration similar to the one 
presented bv lSantos et~al] (|2000h ; 

tt From HA RPS spectra using a calibration similar to the one 
presented bv ISantos et al] (|2002T ) 



The effective temperature and surface gravity derived 
are compatible with the expected parameters of a slightly 
evolved solar- type star, being roughly compatible with 
the spectral type of GO listed in the Simbad dat abase. 
Using the metho d described in Pont fc Ever! (|2004l ) with 
the isochrones of Girardi et al. ( 2000[ ) we derived a stellar 
mass and age of 0.99i0.08M Q and 6-11 Gyr, respectively. 
The stellar radius inferred is 1.95±O.26R0. The uncertain- 
ties in these parameters were derived using 2-sigma errors 
in the effective temperature, surface gravity and metallic- 
ity. Using this mass and radius together with the derived 
surface gravity and effective te mperature, an d making use 



1996), we derive a 



of the bolometric correction of lFlower 
distance of 90 pc to HD 171028, with an error around 15 pc 
if we consider a conservative error in the logg of 0.15 dex. 
No reddening corrections were taken into account in this 
estimate. 

The values menti oned above are a bit different from 
those listed in the iNordstrom et al 
(T eff = 5432 if, [Fe/H] = -0.81 ± 0.02 



(|2004f ) catalogue 
Mass=0.78M Q , 



and distance of 43 pc). We note, however, that the pho- 
tometric calibrations used by Nordstrom et al. to derive 
these parameters (including the absolute magnitude used 
to derive the distance) are likely only valid for main- 
sequence stars. 

From the HARPS spectra we derived a chromospheric 
activity index (logi?y K = —4.92, with rms of 0.02) 
following a procedu re similar to the one presented in 
(jSantos et all l2000h . From the activity level and the 
B — V colour we estimate a rotationa l period of 19 days 
(|Noves et al.lll984l) and an age of 4 Gyr dHenrv et al.l ll996l) 
(or at least above 2 Gyr - IPace fc Pasauini 2004 ) . These 
latter values are likely not accurate due to the fact that 
HD 171028 is a metal-poor s tar and a bit evolved off the 
main sequence ( WrightJ[200J) . The calibrations mentioned 



Table 2. Chemical abundances for the several alpha and 
iron peek elements studied. For titanium, both results 
based on Ti I and Ti II lines are presented. The abundances 
given by Fel and Fen lines is the same (it is one of the 
conditions for the d erivation of the stellar parameters - 
Santos et al. 2004b] ). The error in [X/H] represents the 



rms of the abundances given by the n(X) different lines 
used. 



Element 


loge© 


[X/H] 


n(X) 


[X/Fc] 


Fe 


7.47 


-0.49T0.02 


38 


0.00 


Si 


7.55 


-0.43T0.02 


11 


0.06 


Ca 


6.36 


-0.41i0.05 


13 


0.08 


Sc 


3.10 


-0.40i0.05 


6 


0.09 


Til 


4.99 


-0.44i0.03 


12 


0.05 


Tin 


4.99 


-0.43i0.04 


4 


0.06 


V 


4.00 


-0.54i0.04 


8 


-0.05 


Cr 


5.67 


-0.56i0.02 


4 


-0.07 


Mn 


5.39 


-0.72i0.06 


4 


-0.23 


Co 


4.92 


-0.50i0.03 


6 


-0.01 


Ni 


6.25 


-0.53i0.03 


29 


-0.04 


Na 


6.33 


-0.42i0.07 


3 


0.07 


Mg 


7.58 


-0.44i0.03 


2 


0.05 


Al 


6.47 


-0.46i0.01 


2 


0.03 



tivity level is compatible with HD 171028 being an old, 
chromospherically quiet star. 

3.2.1. Abundances of different species 

It is known that thick-disk stars have typically higher 



abundances of alpha elements (jBensbv et al] 12003 



Fuhrmannl 120041 ). To try to access if HD 171028 could 
be a member of this galactic po pulation, we used the 
method and line-lists described in Gilli et al. (120061 ) and 
Santos et al. ( 20061 ) to derive the abundances of several 
alpha and iron-peek elements for HD 171028. The results, 
listed in Tabled suggest that this star has ty pical abun- 
danc es of thin disk solar-type stars (see also iGilli et all 
2006). In particular, no overabundance of alpha elements 



is seen when comparing with stars of similar [Fe/H]. 

No trace of the Li line at 6707.8 A is found in our 
S/N~500 spectrum of HD 171028. An upper limit of 
0.6 mA was obtained for the Equivalent Width of the 
Li-line. This value translates into an upper limit of 
loge(Li) < 0.2 dex. Such a value is compatible with its 
evolutionary status and effective temperature, although 
a large dispersion is seen in the Li abundances of sub- 
giant stars simil ar to HD 171028 (jRandich et al. 
iLebre et al1ll999h . 



1999; 



above were based on main-sequence stars, mostly of solar 
metallicity. In any case, the measured chromospheric ac- 



3.3. HARPS orbital solution 

As mentioned above, we obtained 19 accurate radial- 
velocity measurements of HD 171028 with the HARPS 
spectrograph, between October 2004 and April 2007. The 
complete radial velocity measurements obtained and the 
corresponding errors are presented in Table[3] It is worth 
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Fig. 2. Top: Phase-folded radial- velocity measurements 
of HD 171028, and the best Keplerian fit to the data 
with a period of 538-days, eccentricity of 0.61, and semi- 
amplitude of 58ms -1 . 

noticing that the errors quoted in the table, which are 
used to plot the error bars, refer solely to the instrumen- 
tal (calibration) and photon-noise error share of the total 
error budget (e.g. activity and/or stellar oscillations are 
not considered, given the difficulty in having a clear esti- 
mate of their influence). 

Just after the first measurements were done the star 
was noticed to be radial-velocity variable. A later analysis 

Table 3. HARPS radial-velocity measurements of 
HD 171028. 



JD 


V T [kms -i ] 


a(V r ) [kms -1 ] 


2453310.4946 


13.6319 


0.0009 


2453574.6582 


13.5948 


0.0008 


2453575.6376 


13.5924 


0.0007 


2453576.6401 


13.5899 


0.0008 


2453577.6084 


13.5939 


0.0008 


2453578.6591 


13.5954 


0.0007 


2453669.5027 


13.7099 


0.0016 


2453672.4942 


13.7045 


0.0013 


2453862.8107 


13.6268 


0.0011 


2453864.8024 


13.6304 


0.0011 


2453870.7806 


13.6268 


0.0009 


2453882.8334 


13.6266 


0.0010 


2453883.8141 


13.6219 


0.0017 


2453920.7281 


13.6179 


0.0008 


2454166.9019 


13.6329 


0.0009 


2454168.8885 


13.6345 


0.0008 


2454172.8740 


13.6417 


0.0008 


2454174.8916 


13.6446 


0.0008 


2454196.8946 


13.7075 


0.0007 



Fig. 3. Top: Radial-velocity measurements of HD 171028 
as a function of time, and the best Keplerian fit to the 
data. Bottom: Residuals of the fit. 



of the whole data set revealed the presence of a 538-day pe- 
riod radial-velocity signal. This signal is best fitted using 
a Keplerian fit with an semi- amplitude K of 58ms -1 , and 
an eccentricity of 0.61 (see Figs.[5]and[3]). Given the mass 
for HD 171028, this corresponds to the expected signal in- 
duced by the presence of a 1.83 Jupiter-masses (minimum- 
mass) companion (Table|4]). 

To understand if the periodic radial-velocity signal 
observed could have a non-planetary or i gin (see e.g. 
Saar fc Donahue! Il997l : iQueloz et all l2000t ISantos et all 
2002), in Fig.g]we plot the Bisector Inverse Slope (BIS) of 
the Cross- Correlation Function as a function of the radial- 
velocity. The result shows that no correlation exist be- 
tween the two quantities, suggesting that stellar activity 
of stellar blends cannot explain the radial-velocity varia- 
tion observed. Together with the low activity level of the 
star, we conclude that the 538-day orbital period observed 



Table 4. Elements of the fitted orbit for HD 171028b. 



p 


538±2 


[d] 


T 


2453648.9205±1.8584 


[d] 


a 


1.29 


[AU] 


e 


0.61±0.01 




V r 


13.631±0.001 


[kms -1 ] 


■jj 


305±1 


[degr] 
[ins -1 ] 


Ki 


58.0±0.4 


Mm) 


5.421 10~ 9 


[M Q ] 


a(0-C) 


2.4 


[ms- 1 ] 


N 


19 




77i2 sini 


1.83 


[M Jup ] 
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RV [m/s] 

Fig. 4. BIS vs. radial-velocity for HD 171028. To strongly 
evidence the nonexistence of a correlation the vertical and 
horizontal scales were set to be the same. 



can be better explained by the presence of a Jupiter-like 
planet orbiting HD 171028. 

The residuals of the orbital fit have a rms of 2.4 ms -1 , 
slightly above the average photon-noise error of the mea- 
surements (0.95 ms -1 ). The lower panel of the plot pre- 
sented in Fig.[3]does reveal some structure in the residuals 
after the 538-day period orbit is subtracted. Although we 
cannot exclude that this signal may be due to the presence 
of another companion to the system, the fact that some 
structure is also present in the plot of Fig. Q] may hint at a 
non-planetary origin for the signal. More data are needed 
to settle this issue. 



4. Concluding remarks 

In this paper we present the detection of the first planet 
among the stars in the HARPS metal-poor planet search 
program, a giant planet with a minimum mass of 1.8Mj up 
orbiting HD 171028 in an eccentric trajectory every 538 
days. 

This detection adds to the small number of planets 
kno wn to orbit stars w i th metallicity clearly below solar 



(e.g. iMavor et al.ll2004 ICochran et alj 120071) . While the 



correlation between the presence o f planets and stellar 



metallicity is clearly established (e.g. [G onzale z et al. 2001 



Santos et al.||200lll2004ri IFischer fe Valentill2005l )~the de 



tection of an increasing number of giant planets orbiting 
low-metallicity stars reopens the debate about the origin 
of these worlds. These findings show that giant planet for- 
mation is not completely inhibited around stars in the 
intermediate metal-poor regime. 



These cases either represent the metal-poor tail of the 
giant planets formed by the core-accretion process, or they 
may hint at the existence of a different population of 
planets formed as a result of disk-instability processes. 
The precise study of the metallicity distribution of stars 
with planets suggest s that there may ex i st a "long" flat 
low-m etallicity tail ( Santos et al.1 l2004bl Udrv fe Santosl 
20071 ). Although the statistics of metal-poor planet-hosts 



is still poor, we can speculate that we may be observ- 
ing a superposition of two populations. On the one hand, 
stars whose planets were for med by the metallicity depen- 
dent c ore-accretion process ([Pollack et al ][l996HIda fc Li3 
2004b), mostly populating the metal-rich regime. On the 



other hand, a less significant population of giant planets 
orbiting stars of all metallicities. These latte r could hav e 
been formed by the disk-instability process (|Bossll2002h . 
We are considering here that subsequent planet evolution 
processes (e.g. migration in the disk) are not stro ngly de- 
pendent on the metallicity (|Livio fc Pringldl2003h . 

We should add, however, that according to the models, 
giant planets formed by the disk-instability process could 
have higher masses w hen compared with those formed 
by core-accretion (e.g. Matsuo et al.l l2007h . The lack of 
a clear co rrelation between p lanet-mass and stellar metal- 
licity (e.g. Santos et al.ll2003h may be an important caveat 
for the proposed scenario. 

In this context it is interesting to verify that HD 171028 
is clearly one of the most metal-rich stars in the HARPS 
sample presented in this paper. Although statistically not 
relevant (note that our sample is clearly more populated at 
the metallicity range between —0.5 and —0.8), if confirmed 
this fact could lend some support to the core-accretion 
model. 

Also interesting is the fact that HD 171028 is a bit 
evolved off the main -sequence. In a very recent study, 
Pasauini et al. (2007) present evidence that giant stars 



with planets are likely not as metal-rich as their dwarf 
counterparts (see also Ida Silva et al. 2006a ). Although 
other explanations exist (e.g. the higher mass of the 
stars, and eventually of the proto-planetary disks, may 
significantly change the planet formation efficiency - 
Laughlin et al _J|2004 lEndl et alJliool iBonfils etai]|2007 : 



Johnson et al.ll200*^ ). this interesting conclusion could sug- 



gest that planetary pollution may be more important 



than previously thought ( e.g. |P insonneault et al] 12001 



Santos et all I2003L l2004ri iFischer fc Valentil 120051) . We 
note, however, that no clear connection seems to be 
present between the stellar evolution status and [Fe/H] 
among dwarfs and sub-giant stars with planets (e.g. 
Santos et all 120031 : IFischer fc Valentil l2005h . In the sub- 



giant branch, planet-host stars are still metal-rich when 
compared with "single" field stars, even though impor- 
tant dilutio n processes may have occu rred. A look at 
the ta bles of lSantos etHl (|2004bL l2005h and ISousa etal 



( 20061 ) , also shows that low metallicity among planet hosts 
is not restricted to evolved stars. Such an observations 
could be expected if stellar pollution was a frequent out- 
come of the planet formation process. Out of 14 stars 
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with metallicity below — 0.20dex, only 5 are likely evolved 
(log<7 <4.2), while 7 "definite" dwarf stars (log g >4.4) 
are also metal-poor. On the other hand, several evolved 
stars exist with metallicity well above so lar. Finally, a l- 



though some caveats have been discussed (|Vauclairll2004l ). 
the lack of any correlation between stellar metallicity 
and convective envelope mass is still an important argu- 
ment against pollution being the main mechanism respon- 
sible for t he observed metal rich n a ture of planet-host 
stars (e.g. IPinsomreault et al. 2001 ; Santos et al.l [2003; 
Fischer fe Valentill2005ir The K-dwarf stars with planets, 
with deep convective regions, would require an excessive 
infall of planetary material to get enriched to the observed 
level. 

The continuation of the program presented here will 
certainly provide important constraints for this discussion. 
On the one hand, it will help to understand what is the 
lower stellar metallicity limit for which giant planet s are 
still able to form (see discussion in lMatsuo et al.ll2007l ). On 
the other hand, with the adequate observing strattegy, it 
will allow to access the frequency of N eptune-mass plan - 
ets orbiting lower metallicity stars (see lUdry et al.ll2006l ). 
Although the (unknown) disk masses may also have a cru- 
cial influence on the planet formation efficiency, such re- 
sults would have clearly important implications for the 
models of planet formation and evolution. 
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